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ABSTRACT 

The success of any ALMA phase-calibration strategy, which incorporates phase 
transfer, depends on a good understanding of how the atmospheric path delay 
changes with frequency (e.g. Holdaway & Pardo 2001). We explore how the wet dis- 
persive path delay varies for realistic atmospheric conditions at the ALMA site using 
the ATM transmission code. We find the wet dispersive path delay becomes a signifi- 
cant fraction (> 5 per cent) of the non-dispersive delay for the high-frequency ALMA 
bands (> 160 GHz, Bands 5 to 10). Additionally, the variation in dispersive path delay 
across ALMA's 4-GHz contiguous bandwidth is not significant except in Bands 9 and 
10. The ratio of dispersive path delay to total column of water vapour does not vary 
significantly for typical amounts of water vapour, water vapour scale heights and 
ground pressures above Chajnantor. However, the temperature profile and particu- 
larly the ground-level temperature are more important. Given the likely constraints 
from ALMA's ancillary calibration devices, the uncertainty on the dispersive-path 
scaling will be around 2 per cent in the worst case and should contribute about 1 per 
cent overall to the wet path fluctuations at the highest frequencies. 



1 INTRODUCTION 

The performance of interferometers at (sub)millimetre wave- 
lengths is often limited by differential fluctuations in the at- 
mospheric path along the line of sight to each of the con- 
stituent antennas. If uncorrected, these fluctuations lead to 
a loss in sensitivity, imaging artefacts and a limit on the max- 
imum usable baseline (e.g. Carilli & Holdaway 1999; Asaki 
et al. 2005; Nikolic, Richer & Hills 2008). ALMA will correct 
such path variations using a combination of techniques: 

(i) Fast-switching observations of bright calibrator sources 
(e.g. quasars). 

(ii) Water vapour radiometry using dedicated 183-GHz ra- 
diometers (WVRs), installed in every 12-m ALMA antenna. 
The WVRs will allow the retrieval of the amount of water 
vapour along the line of sight to each antenna and thus fluc- 
tuations in the atmospheric path length resulting from this 
water vapour. 

(iii) Self calibration, for a fraction of the brightest science 
targets. 

Fast switching interleaves science observations with short (~ 
a few seconds) calibration observations on a cycle time of 20- 
200 s. These calibrations allow an estimate of the atmospheric 



and instrumental phase errors in the direction of the calibra- 
tor. The applicability of the calibrator phase solutions to the 
science target depends on their angular separation and the 
observing frequency. The radiometric technique will be ap- 
plied continuously, correcting for atmospheric path fluctua- 
tions on timescales shorter than the fast-switching cycle time, 
but does require an accurate model to convert the WVR mea- 
surements into path delays. 

Variations in the atmospheric path delay to the anten- 
nas predominantly arise from a combination of fluctuations 
in the water vapour content and density of air in the tropo- 
sphere - respectively the so-called wet and dry path com- 
ponents. Furthermore, we can split both components (the 
wet and dry) into two parts: one dispersive (i.e. dependent 
on frequency) and the other non-dispersive (independent of 
frequency). Conventionally, the non-dispersive part is taken 
to be the path in the low-frequency limit. This memo fo- 
cuses on the dispersive part of the wet path delay intro- 
duced by the atmosphere. Previously atmospheric disper- 
sion in the context of fast switching was studied by Hold- 
away & Pardo (2001) and Holdaway & D'Addario (2004). 
They found the magnitude of the dispersive phase will be- 
come non-negligible in ALMA's submillimetre bands. Fur- 
thermore, Holdaway & Pardo (2001) quantified the disper- 
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sive phase delay for typical conditions at the ALMA site. We 
extend their work by investigating how the dispersive path 
delay depends on the variation of the physical parameters of 
the atmosphere and what constraints can be placed on it from 
the proposed ancillary calibration devices on site. 

This memo is structured into six parts. First, we ex- 
plain the planned phase-correction strategy for ALMA and 
where dispersive effects play a role. Second, in § 3, we de- 
scribe how dispersive and non-dispersive path delays arise 
in the atmosphere and how we compute the dispersive path 
delay. § 3 also quantifies the dispersive path delay across 
the full range of ALMA observing frequencies, before § 4 
demonstrates how varying different atmospheric parameters 
changes the magnitude of the dispersive path delay. Finally, 
we look at what constraints can be placed on the dispersive 
path delay with temperature measurements from weather- 
monitoring devices. 



2 PHASE CALIBRATION FOR ALMA 

Before we consider the magnitude of the dispersive path 
delay, we will detail the current phase-correction plans for 
ALMA and how these are affected by atmospheric disper- 
sion. 

Numerous previous memos have discussed fast- 
switching phase correction for ALMA (e.g. Carilli & Hold- 
away 1999; Holdaway 2001; Holdaway & D'Addario 2004). 
Briefly, fast switching removes some of the antenna-based 
phase fluctuations by regularly making short observations of 
a bright point-source calibrator, close to the science target in 
the sky. This difference in direction between the target and 
calibrator means the calibration solutions need to be interpo- 
lated to estimate the science target's phase errors. 

Typically, at ALMA's highest observing frequencies, it 
will be unlikely that the calibrations can be taken at the same 
frequency as the science observations, for the following rea- 
sons: 

(i) Calibrators of the necessary strength at the science fre- 
quency will probably be too far from the science target (al- 
though this depends on the source counts at high frequen- 
cies, which are currently not well known). 

(ii) The combination of high frequencies and long base- 
lines may resolve many potential calibrators. 

(iii) The phase errors may prove to be so large that phase 
wrapping makes reliable solutions difficult. 

To overcome these difficulties, one proposal is for ALMA 
to observe calibration sources at 90 GHz when necessary 
and then scale the corresponding phase solutions to the sci- 
ence frequency (the phase-transfer technique). The frequency 
above which the phase transfer-approach will be necessary 
depends on the array configuration, atmospheric conditions 
and ultimately experience gathered at the site. It is currently 
expected to be the routine mode at wavelengths below 2 mm 
(frequencies > 150 GHz). In the absence of dispersion, the 
required phase scaling is simply the ratio of the calibration 
and science frequencies. However, at high frequencies, as we 
have already mentioned and will detail below, the numerous 
nearby water vapour lines ensure that atmospheric disper- 
sion will need to be taken into account. 

Fluctuations in the atmospheric path on timescales 



shorter than the fast-switching cycle time will be corrected 
radiometrically using WVRs installed in every 12-m ALMA 
antenna. The WVRs measure the brightness of the 183-GHz 
atmospheric water line, which is very sensitive to the total 
amount of water vapour in the atmosphere (Stirling et al. 
2004; Nikolic 2009a). The fundamental difficulty in the anal- 
ysis of WVR data is how to convert fluctuations in the mea- 
sured sky brightnesses around the 183-GHz water line into 
variations in the atmospheric path delay. We have begun 
to develop a Bayesian framework for computing such con- 
version factors (Nikolic 2009a,b), which naturally incorpo- 
rates the prior knowledge of the system and all the observ- 
able data. In the first memo (Nikolic 2009a), we begin with 
the simplest possible model atmosphere, comprising a sin- 
gle, thin layer of water, which is the only cause of path 
fluctuations. Even so, an analysis of test data from proto- 
type WVRs that were installed on the Submillimeter Array 
(SMA) yields corrections that are within ~ 5 per cent of op- 
timal ones. Memo 588 (Nikolic 2009b) extends this scheme 
to include the observed correlation between phase and sky 
brightness, which ALMA may implicitly record during fast 
switching, when the calibrator's phase is measured whilst 
data are taken with the WVRs. The inclusion of this empir- 
ical relationship should significantly improve the accuracy 
of the phase correction. 

Typically, the calculation of phase-correction coefficients 
(i.e. the conversion factors above) from WVR data requires 
a model of the atmosphere, which in principle can be used 
to estimate the dispersive as well as non-dispersive path de- 
lays. In our work to date, which was based on relatively low- 
frequency (220 GHz) data, we have not modelled the dis- 
persive path contribution. Instead, we simply scale the non- 
dispersive path delay by a constant factor to estimate the 
total wet delay. However, in reality this factor is a function 
of frequency. Furthermore, if the dispersive contribution is a 
significant fraction of the total path delay, it must be deter- 
mined with good accuracy. To date there has been little work 
on what atmospheric properties influence the magnitude of 
the dispersive path delay and by what amount. In § 4, we 
therefore examine the variation of the dispersive path delay 
using multiple realistic models of the atmosphere above the 
ALMA site. 



3 THE DISPERSIVE PATH DELAY 

Fluctuations in the wet path delay, d/ H2 o, are separated into 
the sum of non-dispersive (ds) and dispersive delays (dS v ): 



d'H 2 = ds + dS v . 



(1) 



The total non-dispersive path delay (i.e. the sum of the wet 
and dry components) can be computed from the Smith- 
Weintraub equation for the refractive index, n, at a temper- 
ature, T (Stirling et al. 2008; Nikolic 2009a): 



10" 



4+, 



j + 7 j2 



(2) 



where and P w are the partial pressures of the dry air and 
water vapour respectively and a, fi and 7 are constants. Since 
we are only concerned with the excess path introduced by 
water vapour (i.e. ds above), then we may omit the first term 
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and of the remaining terms the last one dominates. Hence we 
can transform to the following expression (Nikolic 2009a): 



where c is the water vapour column. 

In general, atmospheres with spectral lines are disper- 
sive, with the dispersion related to the absorption accord- 
ing to the Kramers-Kronig relations. Thus, any atmospheric 
property that affects the shapes of absorption-line features 
will probably influence the atmospheric dispersion, summa- 
rized in the following functional form (see Tab. 1 for a list of 
the symbols used): 

d/ H2 o« l^-^)dc + dS v (v,T,T,p,c), (4) 

where we have substituted the non-dispersive delay from 
Eq. 3 into Eq. 1. The variation of dS v as a function of its pa- 
rameters is studied in the remainder of this memo. 

3.1 Computing the dispersive path delay 

We calculate the dispersive and non-dispersive contributions 
to the wet and dry path delays between 1 GHz and 1 THz 
from model atmospheres using the ATM software (Pardo, 
Cernicharo & Serabyn 2001) 1 . ATM accurately predicts the at- 
mospheric opacity above Mauna Kea, Hawai'i up to 1.6 THz 
(Pardo et al. 2005). Its dispersive calculations have not been 
as thoroughly verified due to a lack of high-frequency test 
data, which should change in the near future once the ALMA 
WVRs become operational at Chajnantor. 

In Fig. 1, we plot the wet and dry, dispersive and non- 
dispersive path delays overlaid on top of the atmospheric 
transmission, all computed using the basic model parame- 
ters listed in Tab. 1, which are suitable for the ALMA array 
operations site (AOS) at Chajnantor. The parameters listed 
in Tab. 1 are the median values measured by site-monitoring 
equipment where possible and their origin will be detailed in 
later sections. The computed dispersive path delay is small 
except around the atmospheric H2O lines (for the wet com- 
ponent) or O2 lines (for the dry). The largest contribution to 
the total delay at all frequencies arises from the dry non- 
dispersive path delay. However, we expect the dry column 
to be relatively stable and so contribute very little to the dif- 
ferential path delay (Holdaway & Pardo 2001; Holdaway & 
D'Addario 2004) i.e. the difference in delay between anten- 
nas. We will ignore the dry path delays for the rest of this 
memo but they may prove to be a significant extra source of 
error that remain after phase correction using the WVRs. 

In Fig. 2, we plot the ratio of wet dispersive to non- 
dispersive path delay across the (currently-funded) ALMA 
observing bands (see Tab. 2). In the first two of these bands 
(Bands 3 and 4), the dispersive path delay is between 0.5 and 
3 per cent of the non-dispersive, rising as a fraction approx- 
imately linearly with frequency. At the frequencies of Band 
5 and higher, the dispersive path becomes a more significant 

1 We have packaged ATM as AATM, which is available under 
the GPL license from http://www.mrao.cam.ac.uk/~bn204/alma/ 
atmomodel.html. We interface with the ATM libraries using the 
dispersive program (Nikolic 2009a), from version 0.1 of AATM. 



Table 1. Basic parameters of the atmospheric model of the ALMA 
AOS site, computed from median site characterization data except 
Pg (see Appendix A for details). 



Parameter 


Units 


Value 


Comment 


c 


mm 


1.22 


Zenith water column 


Pg 


mb 


560 


Ground-level pressure 


T , 


K 


270 


Ground-level temperature 


r T 


Kkm- 1 


-7.28 


Tropospheric lapse rate 


ho 


km 


1.16 


Water vapour scale height 



Table 2. ALMA receiver bands. In the initial phase of operations 
ALMA will be equipped with only six bands: 3,4, & 6-9. Bands 1 
and 2 are not yet funded. 



Band 


Frequency 


Trx 


v 1 
u rep 


Mixing 3 


Supplier 4 




range (GHz) 


(K) 


(GHz) 


Scheme 




1 


31-45 


17 


38 


USB 




2 


67-90 


30 


79 


LSB 




3 


84-116 


37 


100 


2SB 


HIA 


4 


125-163 


51 


144 


2SB 


NAOI 


5 


163-211 


65 


200 


2SB 


osot 


6 


211-275 


83 


243 


2SB 


NRAO 


7 


275-373 


147 


342 


2SB 


IRAM 


8 


385-500 


196 


405 


2SB 


NAOI 


9 


602-720 


175 


680 


DSB 


NOVA 


10 


787-950 


230 


869 


DSB 


NAOI 



1 Receiver noise temperature specification for over 80 per cent of the 
band. 



2 Representative frequency, either the band centre or where the trans- 
mission is better if the centre is near an absorption line. 

3 Two lowest-frequency bands use HEMT mixer technology and are 
single sideband, either upper (USB) or lower (LSB); all the others 
use SIS mixers and are either dual sideband (2SB - each sideband 
detected separately) or double sideband (DSB). 

4 HIA - Herzberg Institute of Astrophysics, Canada; NAOI _ Na- 
tional Astronomical Observatory of lapan; OSO - Onsala Space Ob- 
servatory/ Chalmers University, Sweden; NRAO - National Radio 
Astronomy Observatory, USA; IRAM - Institut de Radio Astronomie 
Millimetrique, France; NOVA - Nederlandse Onderzoekschool voor 
de Astronomie. 

+ Six production receivers will be provided through the European 
Commission's Framework 6 programme. 

fraction of the total wet delay (> 5 per cent) and will need to 
be considered in both the fast-switching and WVR analyses. 
In Band 8 for instance, as far as possible from the absorption 
lines, the dispersive path delay is 20 to 40 per cent of the non- 
dispersive. 

Additionally, we can investigate the phase slope that 
the atmosphere would introduce across the observing band- 
pass if no account were taken for variations in the dispersive 
path delay with frequency. In Fig. 3, we plot the fractional 
variation in the dispersive path delay for a 4-GHz portion 
of frequency space around representative band frequencies 
listed in Tab. 2. For each band we have normalized the val- 
ues with respect to the chosen bandpass centre frequency. 
As we noted for Fig. 2, the dispersive path delay varies ap- 
proximately linearly for small frequency widths. For most of 
the bands the variation is ±2 — 5 per cent over the observ- 
ing bandwidth. Thus, such variation can be ignored, particu- 
larly at the lowest frequencies where the dispersive path de- 
lay makes a small contribution overall. Of the bands plotted, 
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Figure 1. Wet non-dispersive (fop panel, dashed blue) and dispersive (top panel, solid blue) path delays alongside the dry non-dispersive (middle 
panel, dashed red) and dispersive (bottom panel, solid red) path delays calculated by ATM for the AOS model atmosphere with standard param- 
eters as listed in Tab. 1, including a 1.22 mm column of water. Overlaid in grey on all three panels is the corresponding zenith atmospheric 
transmission, calculated from the total opacity x (as e~ T ). 
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Figure 2. Ratio of the wet dispersive to wet non-dispersive delay for the construction ALMA observing bands, listed in Tab. 2. Also plotted {grey) 
is the atmospheric zenith transmission at the same frequency. The ratio rises continually with frequency and becomes particularly significant in 
the submillimetre bands (Bands 7,8 & 9). 



only Band 9 exhibits large variations of some ±13 per cent 
across the 4-GHz chunk, which rises to nearly ±30 per cent 
if the full 8-GHz DSB bandwidth is considered. In this case, 
a frequency-dependent scaling or multiple independent cal- 
culations of the dispersive phase (as planned in the ALMA 
correlator and TelCAL software sub-system) will be impor- 
tant. 



4 PHYSICAL INFLUENCES ON THE DISPERSIVE 
PATH DELAY 

In this section we quantify the impact of changes in atmo- 
spheric conditions on the wet dispersive path delay. We ex- 
plore the influence of: the quantity of water vapour (§ 4.1), 
the ground-level temperature (§ 4.2) and pressure (§ 4.3) 
alongside the distribution of temperature with height (§ 4.4) 
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Figure 3. Variation in wet dispersive path delay for the construction 
ALMA bands, which were listed in Tab. 2. The lines use the same 
colours as Fig. 2 for the different bands: 3 (purple), 4 (dark blue), 5 
(light blue), 6 (green), 7 (yellow), 8 (orange), 9 (red). The delays are plot- 
ted for frequencies in a 4-GHz portion around the bands' represen- 
tative frequencies (listed in Tab. 2, typically the band centres). This 
portion represents half of the instantaneous contiguous bandwidth 
(either the upper or lower sideband) for the ALMA 2SB mixers. Each 
measurement is divided by the value at the representative frequency 
and all the computations are for the median atmosphere specified in 
Tab.l. The Band 9 receivers are DSB and therefore provide 8-GHz 
contiguous bandwidth (not plotted), which results in about a 30 per 
cent variation in the dispersive path delay towards the band edges 
compared with the representative frequency. 

and the distribution of water vapour with height (§ 4.5). For 
each investigation we compute the dispersive delay from 
ATM using a variety of atmospheric models designed to rep- 
resent average and extreme conditions at the ALMA AOS. 

Our focus is the radiometric phase-correction technique 
using the WVRs. An important requirement for this tech- 
nique is the optimal computation of time-dependent phase- 
correction coefficients, dL/dT B ,', for each of the four WVR 
channels. These coefficients relate a change in atmospheric 
path, SL, to a change in the observed sky brightness by a 
WVR, in each of its four channels, <5T B via (see Nikolic 
2009a,b): 

4 dL 

SL='£zv l —ST B „ (5) 

i=l dT B,< 

where Wj is an appropriate weighting for each channel. We 
can decompose this expression further, since the WVRs are 
really only sensitive to variations in the water vapour col- 
umn, 8c, which we can relate to the measured sky brightness 
by introducing coefficients, Ac I dT B/ ;: 

4 Ac 

Sc=Y^Wi--—ST B:l . (6) 

Although the WVRs do place constraints on the atmospheric 
temperature and pressure, normally they cannot detect small 
fluctuations in their values. We can then write the variation 
in wet path delay, SI H2 q, as the product of the fluctuation in 
water vapour content and the sum of two scaling terms: 

/1741K dS v \ 
Sl H20 =^^^ + —JSc. (7) 

The first scaling term is the non-dispersive path delay from 



Eq. 3, while the second quantifies the dispersive path delay. 
Finally, as we saw in § 3, although AS V / Ac may vary a lot, it is 
only related to the dispersive path delay, which in the millime- 
tre bands makes a very small contribution to the total path 
delay compared to the non-dispersive path. 

4.1 Water vapour quantity 

First, we investigate the effect of varying the column of water 
vapour, c, in the standard AOS atmosphere (Tab. 1). In Fig. 4 
we plot AStVAc for c = 0.44, 0.69, 1.22, 2.56 and 5.45 mm cor- 
responding to the 10, 25, 50, 75 and 90 percentiles of the pre- 
cipitable water vapour (PWV) cumulative function at Cha- 
jnantor respectively. These percentiles are derived from the 
cumulative 225 GHz opacity distributions as described in 
§ Al of Appendix A. Fig. 4 is the template for all the plots 
that follow in this section. The upper panels display the ab- 
solute variation in AS V / Ac for the full extent of the ordinate 
(left panel) and for an enlarged portion close to the horizon- 
tal axis (right panel). In the lower panels, we plot the relative 
change in AS V / Ac, computed by dividing its spectrum (i.e. 
AS V I Ac as a function of frequency) by the one from the me- 
dian column, c = 1.22 mm. Again, on the left-hand side we 
show the full ordinate range, whilst on the right we enlarge 
the plot around a ratio of unity. 

The shape of the AS v /Ac spectrum closely resembles 
that of the dispersive path delay overall, i.e. AS V is propor- 
tional to Ac (see Fig. 1), having discontinuities where the dis- 
persive path delay wraps around, near the dips of the atmo- 
spheric transmission and bright water lines. Even in the en- 
larged (upper right panel) and relative plots (lower panels) there 
are no discernible differences between AS V / Ac for the differ- 
ent values of c. The only significant departures from unity for 
the relative AS v /Ac are at discontinuities in AS v /Ac. Cor- 
respondingly over the entire range of conditions likely at 
the ALMA AOS, the amount of water vapour does not af- 
fect ASv/Ac, only supplying a non-dispersive linear change. 
Furthermore, varying the telescope's elevation /line of sight 
airmass has an equivalent effect to altering the water vapour 
column at zenith. Therefore our plots also indicate that the 
dispersive phase scaling does not depend on airmass and 
will not need to be recomputed for changes in elevation dur- 
ing observing. 

4.2 Air temperature at ground level 

Next, we look at how the ground-level temperature affects 
the predicted phase-correction coefficients. As in the previ- 
ous section (§ 4.1), we varied T g , in the standard AOS atmo- 
sphere which was used as the input to ATM. Fig. 5 provides 
our results for AS V / Ac, in the same format as Fig. 4, using 
T g = 262, 265, 270, 276 and 281 K corresponding to the 10, 
25, 50, 75 and 90 percentiles of the temperature cumulative 
function at Chajnantor (Radford 2004). 

The key plot is in the lower right panel, where we 
plot the spectrum for each T g relative to 270 K. The relative 
AS V I Ac computed are constant with frequency except at the 
centres of absorption features, i.e. the same frequencies as in 
Fig. 4. The variation in the ratio at different Tg is higher than 
for the different water vapour columns, and therefore will 
have more significance. Between the 25 and 75 percentiles, 
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Figure 4. Variation of ASy/Ac with column of water, c. The upper panels show the absolute variation in AS„/Ac as a function of frequency 
from 1 to 1000 GHz, while the lower panels show the variation relative to (i.e. divided by) the median profile with c = 1.22 mm. The right-hand 
panels show the same data as in left-hand ones but over a narrower range of values. The plotted data come from profiles containing 0.44, 0.69, 
1.22, 2.56 and 5.45 mm of water respectively corresponding to the 10, 25, 50, 75 and 90 percentiles of the PWV cumulative function at Chajnantor 
(see § Al). All the curves are plotted in each of the panels: if only one line is discernible then the curves overlap. 



the difference in AS v /Ac is 5-6 per cent but at the 10 and 
90 percentiles it becomes 7-9 per cent. These differences are 
likely to be non-negligible for phase correction at submil- 
limetre frequencies. 

This variation is also important if we recall the strong 
diurnal change in temperature at Chajnantor. Stirling et al. 
(2006) found the temperature varied between approximately 
—5 and 5°C at their site A over 24 hours. Such a change could 
alter the dispersive path delay scaling by around ±5 per cent. 
We return to the constraints we can place on this source of 
uncertainty using meteorological data from the ALMA ancil- 
lary calibration devices in § 5. 



is through the tropospheric lapse rate, Tj. We inves- 
tigate the effect of variations in T T on the dispersive 
path delay in Fig. 7, where AS v /Ac is plotted for Tj = 
-4.80, -5.69, -7.28, -8.83, -9.71 K km" 1 . These Tj again 
correspond to the 10, 25, 50, 75 and 90 percentiles of the T T 
distribution, which we derive from radiosonde data as ex- 
plained in § A2. Realistic variations in the lapse rate at the 
AOS will produce small, non-negligible changes in the dis- 
persive path scaling of some 2-3 per cent from median val- 
ues. Later in § 5, we look at how well we can constrain the 
lapse rate and thus the path scaling using the proposed at- 
mospheric temperature profiler for ALMA. 



4.3 Air pressure at ground level 

We have checked the impact of changing the air pressure at 
ground level, p g , by repeating the previous ATM calculations 
using our standard AOS atmosphere with p g set to 520, 540, 
560, 580, 600 mb. The results, shown in Fig. 6, indicate that 
ASy/Ac varies p s at a level typically less than 0.2 per cent. 
This is much smaller than the effect resulting from any of 
the other atmospheric parameters we consider, so pg is not 
a likely source of uncertainty in the WVR phase-correction 
coefficients. 

4.4 Tropospheric lapse rate 

The main way in which we parameterize the verti- 
cal temperature distribution of the model atmospheres 



4.5 Scale height of atmospheric water vapour 

The final parameter we interogate is the water vapour scale 
height, h . ASy/Ac is plotted in Fig. 8 for h = 0.97, 1.06, 
1.16, 1.29, 1.54 km corresponding to the 10, 25, 50, 75 and 
90 percentiles of the water vapour scale height cumulative 
function at Chajnantor (see § A2). The relative dependence of 
ASy/Ac on hg is again similar to the other parameters, i.e. 
it has spikes at the centres of absorption features. The rela- 
tive AS v /Ac typically varies from the median spectrum by 
between 0.5 and 1 per cent, which rises to 2-A per cent in the 
10 and 90 percentile cases. Thus, under the typical range of 
conditions at the site, ho is a small source of uncertainty in 
the dispersive scaling for phase correction. 
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Figure 5. Variation of AL/ Ac with ground temperature, T g . Panels are laid out as for Fig. 4. The plotted data have T g of 262, 265, 270, 276, 281 K 
corresponding to the 10, 25, 50, 75 and 90 percentiles of the temperature cumulative function at Chajnantor. The temperature data are from 
Radford (2004). 
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Figure 6. Variation of AL/ Ac with ground pressure, p g . Panels are laid out as for Fig. 4. The plotted data have p g of 520, 540, 560, 580, 600 mb. 
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Figure 7. Variation of &L/ Ac with tropospheric lapse rate, Ft. Panels are laid out as for Fig. 4. The plotted data have Tj of —4.80, —5.69, —7.28, 
—8.83, — 9.71 Kkm -1 corresponding to the 10, 25, 50, 75 and 90 percentiles of the Tj cumulative function at Chajnantor (see § A2). — lO.OKkm -1 
is the dry adiabatic lapse rate found in very dry atmospheres. 
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Figure 8. Variation of AL/Ac with water vapour scale height, hg. Panels are laid out as for Fig. 4. The plotted data have ftp of 0.97, 1.06, 1.16, 1.29, 
1.54 km corresponding to the 10, 25, 50, 75 and 90 percentiles of the water vapour scale height cumulative function at Chajnantor. 
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Figure 9. Constraints placed on the dispersive AS v /Ac from a mea- 
surement of Tg. The lines plot the relative change in AS v /Ac from 
the median T g (270 K) spectrum, when T g is changed by 1 K (green) 
or 2K (red). 
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Figure 10. Constraints placed on the dispersive AS v /Ac from a mea- 
surement of lapse rate by the atmospheric temperature profiler. The 
lines show the relative change in AS v /Ac from the typical Tj = 
— 7.28 Kkm -1 spectrum, when Tj is changed by 0.8 Kkm -1 (green) 
or 1.5Kkm _1 (red). 



5 CONSTRAINTS FROM THE ANCILLARY 
CALIBRATION DEVICES 

Our predictions from the atmospheric models above indicate 
that the natural variation of certain atmospheric properties 
at the ALMA site could introduce non-negligible errors in 
the estimation of phase-correction coefficients and also in the 
phase-transfer scaling. However, ALMA will be equipped 
with a suite of ancillary calibration devices, including five 
meteorological towers and a temperature profiler, which will 
monitor basic atmospheric parameters such as the pressure, 
temperature and wind speed. In this section, we look at what 
constraints the temperature devices should be able to place 
on the dispersive path delay prediction. 

First, we examine the ground temperature measure- 
ment. Although temperature probes will be positioned on 
every meteorological tower and possibly on some of the an- 
tennas, we are unlikely to be able to measure T g very accu- 
rately for each antenna. This is mainly because the temper- 
ature profile for the first 100 m of the atmosphere above the 
ALMA site is strongly controlled by surface heating and cool- 
ing, which results in variations of ±5 K over the site (Stirling 
et al. 2006). Therefore, we estimate we will be able to measure 
Tg to around ±2 K for each antenna. 

Using the standard AOS atmosphere with T g = 270 K 
as the input to ATM, we computed the relative change in 
AS„ / Ac for a 1 and 2 K change in Tg, which we plot in Fig. 9. 
At frequencies where the dispersive phase makes a signifi- 
cant contribution to the overall path delay (> 345 GHz), an 
estimate of the ground temperature to ±2 K can constrain the 
dispersive phase-correction coefficient to ±1.5 — 2.0 per cent. 
This is just under the limit of what would be acceptable in 
a total phase calibration budget of 2 per cent. If we could do 
better and constrain Tg to ±1 K, then we could get to better 
than a 1 per cent error. 

ALMA also plans to utilize an atmospheric temperature 
profiler, which will use multi-frequency observations of the 
sky brightness around the 60 GHz O2 lines to infer the atmo- 
spheric temperature profile. Such a profile will provide in- 
formation about the air temperature away from the ground 
where surface effects are minimized and the values are prob- 
ably more settled. The current specifications of the profiler 
unit state it will measure the atmospheric temperature to 



AT ^ 2K with a vertical resolution of Az ^ 200 m up to 
1500 m above the AOS. We assume the temperature profile 
consists of 8 measurements of the temperature separated ver- 
tically by 200 m, and each accurate to ±2 K. If the profile is a 
straight line, we can perform a least-squares fit, which should 
measure the lapse rate, T, to ±1.5 K km -1 . If the accuracy of 
each temperature measurement is better, say to ±1 K, then 
the corresponding accuracy on T reduces to ±0.8 K. Taking 
these constraints on the lapse rate in the standard AOS at- 
mosphere, i.e. Tj = —7.28 ± 1.5 K km -1 , we plot the relative 
change in AS y /Ac in Fig. 10. The plot looks very similar to 
Fig. 9. The error on AS v /Ac in the submillimetre windows, 
where the dispersion is significant, is between 1 and 2 per 
cent using the information from the temperature profiler. If 
the accuracy of the profiler exceeds specifications to mea- 
sure the temperature to ±1 K then this constraint reduces to 
around 0.5-1 per cent. 



6 SUMMARY 

Both of the primary phase-calibration techniques for ALMA 
rely on accurate estimates of how the atmospheric path fluc- 
tuations vary with frequency, because: 

(i) The phase solutions need to be scaled from the cal- 
ibration to the science frequency for phase-transfer fast- 
switching. 

(ii) The WVRs essentially measure variations in the quan- 
tity of water vapour, which must be converted into fluctua- 
tions in the phase delay of the incoming signal. 

First, we revisited the magnitude of the wet dispersive 
path delay in a median AOS atmosphere using the ATM soft- 
ware: 

• In Bands 3 and 4 (84-163 GHz), the wet dispersive path 
delay is a small fraction, 0.5-3 per cent, of the wet non- 
dispersive path delay. 

• For Bands 5 and above (163-720 GHz), the wet disper- 
sive path delay becomes a significant fraction of the non- 
dispersive, > 5 per cent and should be considered in any 
analysis, particularly at the highest frequencies. In the worst 
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case, Band 8 (385-500 GHz), the dispersive path delay is 20- 
55 per cent of the non-dispersive. 

• The variation in the dispersive path delay across the 4- 
GHz instantaneous bandwidth of a single observation is typ- 
ically 2-5 per cent but does rise to higher fractions, ~ 13 per 
cent, in Band 9. Thus, the capability of the ALMA correlator 
and software to apply phase corrections channel-by-channel 
will prove useful. 

Next, we investigated how the the dispersive path de- 
lay changes when the model parameters were varied, over 
ranges that represented the typical and extreme atmospheric 
conditions measured above the ALMA site: 

• The amount of atmospheric water vapour or equiva- 
lently the airmass does not affect AS„/Ac, i.e. the dispersive 
path contribution to the fluctuations. 

• The typical changes in the water vapour scale height and 
ground pressure do produce small changes (< 2 per cent) to 
AS^/Ac. 

• The dispersive path delay depends more strongly on 
the temperature profile of the atmosphere, particularly the 
air temperature at the ground. The 10 to 90 percentiles of 
the Chajnantor ground-temperature distribution cause vari- 
ations in the dispersive AS„/ Ac of 7-9 per cent from the me- 
dian. Additionally, typical diurnal variations in temperature 
(±5 K) would produce similar changes of ±5 per cent. This 
will be significant at frequencies where the dispersive path 
delay provides a major contribution to the total path delay 
(> 345 GHz). ASi//Ac also depends on Yj with the typical 
variation being some 2-3 per cent. 

These results indicate that obtaining a ground-level air 
temperature estimate for each antenna to an accuracy of 
about ±2K will reduce uncertainties in the models of the 
dispersive phase to a satisfactory level. In combination with 
lapse rate estimates, our calculations suggest that the ancil- 
lary calibration instruments on site should be able to con- 
strain the dispersive terms to 1-2 per cent, on target for the 
current calibration budget. 

The relatively small expected variation of dispersive 
path scaling with the natural range of atmospheric condi- 
tions at the AOS is also encouraging from the perspective 
of empirical models, which use the observed correlations be- 
tween phase and WVR measurements. Because the variation 
is small, it means that the frequency-dependence of such em- 
pirical models will not need to be re-calibrated very often. 
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Table Al. Percentiles of the cumulative T225 distribution above Cha- 
jnantor presented in previous work. 



Percentile 


T225 1 


T225 2 


T225 3 


T225 4 


25 


0.036 


0.037 


0.037 


0.037 


50 


0.061 


0.061 


0.062 


0.060 


75 


0.115 


0.122 


0.125 


0.118 



1 Memo 334 (Radford & Chamberlin 2000), spanning 04/95 - 07/00. 

2 Memo 471 (Evans et al 2003), spanning 04/95 - 08/01. 

3 Memo 512 (Otarola et al. 2005), spanning 04/95 - 08/03. 

4 The latest publicly-available data (Radford 2004), spanning 04/95 
- 12/04. 



Table A2. Percentiles of the cumulative T225 and c distributions 
above Chajnantor used in this memo, c has been calculated from T225 
using Eq. Al. 



Percentile 


T 225 1 


c (mm) 


10 


0.026 


0.44 


25 


0.037 


0.69 


50 


0.060 


1.22 


75 


0.118 


2.56 


90 


0.244 


5.45 



1 Latest publicly-available site characterization data (Radford 2004), 
spanning 04/95 - 12/04. 



ALMA Memo 541, The ALMA Project 
Stirling A., Richer J., Hills R., Lock A., 2008, Turbulence 
simulations of dry and wet phase fluctuations at Chaj- 
nantor. Part I: The daytime convective boundary layer., 
ALMA Memo 517.1, The ALMA Project 



APPENDIX A: DISTRIBUTIONS OF 
SITE-CHARACTERIZATION PARAMETERS 

Al Water vapour column 

The quantity of water vapour above Chajnantor is not a 
directly-observable quantity. Instead, we calculate it from the 
measured atmospheric opacity using an appropriate scaling 
for the site. NRAO has operated an automated tipping ra- 
diometer at 225 GHz on the Chajnantor plateau since 1995 
and the data up to August 2004 are collected online (Rad- 
ford 2004). We summarize the various cumulative distribu- 
tions for T225 used in the ALMA memo series and online in 
Tab.Al. 

Various memos have also computed the conversion re- 
lation between T225 and the PWV column (e.g. Delgado et al. 
1999; Bustos et al. 2000). We use the relation presented in Gio- 
vanelli et al. (2001): 

t 225 = 0.0435(PWV/mm) + 0.0068, (Al) 

which is derived for Chajnantor from a comparison of T225 
to the PWV derived from a 183-GHz WVR and agrees well 
with radiosonde data for PWV|3 mm. Using the standard at- 
mospheric parameters we presented in Tab. 1, we also calcu- 
lated the following conversion from ATM: 

T 225 = 0.0416(PWV/mm) + 0.0120. (A2) 

This is similar to Eq. Al but we prefer the experimentally- 
derived scaling as it has been tested more thoroughly. In 
Tab. A2 we list the percentiles of the PWV distribution which 
we use in this memo, derived from the most recent cumula- 
tive T225 distribution that is publicly available (Radford 2004) 
using Eq. Al. 



Table A3. Percentiles of the cumulative ho distribution used in 
this memo, constructed by Bryan Butler who fitted an expo- 
nential function to the data between and 10 km above Cha- 
jnantor (see http://www.tuc.nrao.edu/alma/site/Chajnantor/ 
instruments/radiosonde). 



Percentile 


h (km) 


10 


0.97 


25 


1.06 


50 


1.16 


75 


1.29 


90 


1.54 



of the library of radiosonde launch data above the Chaj- 
nantor plateau. These radiosonde flights were jointly oper- 
ated by Cornell University, NRAO, ESO and the Smithso- 
nian Astrophysical Observatory between October 1998 and 
December 2001. The data and an analysis of each dataset 
are available on dedicated web pages 2 . A preliminary anal- 
ysis of these data was presented in Giovanelli et al. (2001). 
They show the vertical distribution of water vapour den- 
sity in a median atmosphere derived from 108 launches is 
well approximated by an exponential with scale height hy = 
1.135 km. We present the cumulative distribution of hg in 
Tab. A3. This is constructed using the analyses of Bryan But- 
ler, which fitted an exponential to the water vapour pressure 
between and 10 km above the Chajnantor plateau for each 
radiosonde launch. We include 194 separate launches over 
varying months and times of day. 

Furthermore, to ascertain the distribution of the tropo- 
spheric lapse rate, Yj (Tab. A4), we fit a straight line to the ra- 
diosonde temperature data from 204 launches using a least- 
squares method. As most of the water vapour is concentrated 
in the first layers of the atmosphere we fit only to data in the 
first 1 km above the plateau surface. We note that in reality 
few temperature profiles from radiosonde data are a perfect 
straight line and some show more complicated features such 
as temperature inversions. However, for the purposes of this 
memo a simple fit should be sufficient. 



A2 Parameters of the atmospheric vertical profile 

To derive distributions of atmospheric parameters which de- 
pend on the variation of water vapour pressure and tern- 2 http://www.tuc.nrao.edu/aima/site/Chajnantor/ 
perature with height through the atmosphere, we make use instruments/radiosonde 
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Table A4. Percentiles of the cumulative Tj distribution above 
Chajnantor used in this memo. It is constructed from fits to 
the temperature data ^ 1 km from the Earth's surface at Cha- 
jnantor (see http://www.tuc.nrao.edu/alma/site/Chajnantor/ 
instruments/radiosonde). 

Percentile Tj (Kkm -1 ) 
10 -9.71 



25 
50 
75 
90 



8.83 
7.28 
5.69 
4.80 
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